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Abstract 

Refined  estimates  for  the  errors  in  eigenvalue  and  eigenvec¬ 
tor  approximation  by  finite  element,  or,  more  generally,  Galerkir. 
methods,  as  they  apply  to  self-adjoint  problems,  are  presented. 
Particular  attention  is  given  to  the  case  of  multiple  eigenvalues. 
The  results  are  new  in  this  case.  The  proof  is  based  on  a  novel 
approach  which  yields  the  known  results  for  simple  eigenvalues  in 
a  simple  way.  Numerical  computations  are  presented  and  analyzed 
in  light  of  the  theoretical  results. 


I  Acer  !  or  For 

I  N:\s  mUI 


C  o  J  e,  a 


Introduction . 


1 . 

It  is  the  purpose  of  this  paper  to  derive  some  refined  esti¬ 
mates  for  the  errors  in  eigenvalue  and  eigenvector  approximation 
by  finite  element,  or,  more  generally,  Galerkin  methods,  as  they 
apply  to  self-adjoint  problems.  The  results  are  new  in  the  case 
of  multiple  eigenvalues.  The  proof  is  based  on  a  novel  approach 
which  yields  the  known  results  for  simple  eigenvalues  in  a  simple 
way . 

Suppose  is  an  eigenvalue  of  multiplicity  q  of  a  self- 

adjoint  problem  and  let  M(\.  )  denote  the  space  of  eigenvectors 
corresponding  to  '.  ,  where  il  •  I!  denotes  the  energy  norm  for  our 
problem.  Let  S  be  the  finite  dimensional  approximation  space 
employed  in  the  Galerkin  method.  i ^  will  be  approximated  from 
above  by  q  of  the  Galerkin  approximate  eigenvalue: 

'k  "  XS,k' ' S , k+1 . VS,k+q-l’ 

If  we  choose  the  i ^  .  in  increasing  order  then  we  have 

' k  's,k  ‘  ' ‘ ‘  '  '  S , k+q-1 - 

Our  main  estimate  for  the  error  in  eigenvalue  approximation  is 

inf  inf  l|u-i  II  !  2  =  O  (  S  )  2  , 
u-:M(  i,  )  i-:S  I  k 

I 

il  u  II  —  1  j 

thus  showing  that  the  error  between  ' .  and  ^  .  ,  the  appro::!- 
mate  eigenvalue  closest  to  \  ,  is,  to  within  a  multiplicative 
constant,  the  square  of  the  minimal  energy  norm  distance  between 
exact  eigenvectors  u  -  M('^)  with  ;! u !!  =  1  and  S,  i.e.,  the 
square  of  the  energy  norm  distance  between  S  and  the  eigenvector 


u  £  M('^)  with  }! u !!  =  1  that  can  be  best  approximated  by  S. 

For  ,  ,  ,  -  v  ,  the  error  between  and  i_  ,  , ,  the 

S , k+q- 1  k  k  S,k+q-l 

approximate  eigenvalue  farthest  from  ,  we  prove 

(1.2)  A  -  A  i  cf  sup  inf  II u—  i  II ^  2  =  C~  (S)2 

q  1  k  ueMOk)  teS  k 

.  11*11=1 

and  for  the  errors  A  ^  -  v  ^ ,  i  =  we  obtain  bounds 

2 

in  terms  of  quantities  intermediate  in  size  between  Cf  (S)  and 

k 

C',  (S)2. 

k 

These  results  should  be  contrasted  with  those  in  the  litera¬ 
ture.  In  Babuska  and  Aziz  [1],  Fix  [4],  and  Kolata  [6],  the 
estimates 


1.3) 


'  S,k+i  Vk  1  C,';  'k(S>  '  1  0 . q  1 ' 


are  proved.  For  i  =  0,...,q-2f  (1.3)  is  weaker  than  the  esti¬ 


mates  stated  above  ((1.1)  for  i  =  0  and  those  mentioned  after 


(1.2)  for  i  =  i,...,q-2);  for  i  =  q-1,  (1.3)  is  the  same  as 

(1.2)  In  Birkhoff,  de  Boor,  Swartz,  and  Wendroff  [2],  which  pre¬ 
sents  the  earliest  results  of  the  general  type  we  are  discussing, 
the  eigenvalue  estimates  depend  on  the  sum  of  the  squares  of  the 
energy  norm  distances  between  S  and  the  unit  eigenvectors  asso¬ 
ciated  with  all  the  eigenvalues  not  exceeding  i.  .  The 

J 

feature  of  (1.1)  that  is  new  is  the  dependence  of  f  (S)  on  only 

k 

one  eigenvector  u  £  M('^),  namely  the  one  best  approximated  by  S 
Regarding  the  errors  in  the  approximate  eigenvectors,  we  show 
that  if  u^  j.  is  the  Galerkin  approximate  eigenvector  correspond- 
then  there  is  a  u^  =  u^(S)  --  M(*k)  with  ilu^jl  =  * 


ing  to 


'  S  ,  k  ' 


such  that 


<1.3)  II  u  "  u,J  -  Cf  .  (S)  . 

J  t  »\  *  « 

k 

The  error  II  u  .  ,  -  u.  ,  II  is  bounded  by  Cf  ,  (S)  and  the 

S,k+q-j.  k+q-1 

errors  ||us  k+i  -  uk+iH-  1  =  l,...,q-2,  are  bounded  by  quantities 

intermediate  in  size  between  (S)  and  f  (S).  The  best  pre- 

k  '  k 

viously  known  result  is 


(1.4) 


u 


S  ,  k+i 


u 


k+i 


-  Cf 


(S 


i  =  0 


»q-i 


In  Section  2  we  introduce  the  class  of  variational!’/  formu¬ 
lated,  self-adjoint  eigenvalue  problems  considered  in  the  paper, 
define  the  Galerkin  approximations  to  these  problems,  and  in  Lem¬ 
mas  2.1  -  2.3  give  the  preliminary  results  which  are  used  in  the 
sequel.  The  main  theoretical  result  of  the  paper  is  presented  and 
proved  in  Section  3.  The  treatment  is  direct  and  self-contained, 
relying  on  a  minimal  amount  of  functional  analysis  background.  In 
Section  4  we  present  numerical  computations  for  a  finite  element 
approximation  of  a  problem  with  double  eigenvalues  for  which  each 
double  eigenvalue  has  associated  eigenvectors  of  strikingly 
different  approximation  properties.  The  quantities 

2 

sup  inf  ||  u- t  II 

u-=  M  (  V  k )  i  i  S 

nun -i 

and 

2 

inf  sup  || u—  j  || 
usM(>k)  t^S 

II  ull  =  1 

are  thus  of  different  sizes  and  we  would  therefore  expect  '  ,  -v 

S ,  k  k 

and  'g  k+l~'k  to  different  sizes.  This  is  clearly  shown  by 


the  computations.  The  computations  also  show  that  u  (the 

o  /  iC 

approximate  eigenvector  belonging  to  the  approximate  eigenvalue 
closest  to  \^)  converges  to  an  exact  eigenvector  with  good 
approximation  properties,  while  us  k+1 -  the  approximate  eigen¬ 
vector  belonging  to  the  approximate  eigenvalue  farthest  from 

iC 

converges  to  an  exact  eigenvector  with  poor  approximation 
properties . 

The  literature  on  eigenvalue  problems  is  extensive,  with  ma 
papers  bearing,  at  least  tangentially,  on  the  problem  addressed 
this  paper.  We  have,  however,  mentioned  only  those  papers  that 
bear  directly  on  the  central  theme  of  our  results;  namely,  the 
Galerkin  approximation  of  eigenpairs  corresponding  to  multiple 
eigenvalues.  For  a  general  treatment  of  eigenvalue  problems  and 
their  literature,  we  refer  to  the  excellent  and  comprehensive 
monograph  of  Chatelin  [3]. 


Preliminaries . 


2  . 

Suppose  H  is  a  real  Hilbert  space  with  inner  product  (*,*] 
and  norm  11*11,  respectively,  and  suppose  we  are  given  two  symmet 
ric  bilinear  forms  BQ(u,v)  and  D(u,v)  on  H*H.  BQ(u,v)  is 
assumed  to  satisfy 


(2.1) 

and 

(2.2) 


BQ(u,v)  |  ,  C  x  ||  u  ||  II  v  ||  ,  V  u,v  c  H 


C2  II u | 


BQ(u,u),  V  u  a  H,  with  C2  >  0 


It  follows  from  (2.1)  and  (2.2)  that  (u,v)n  =  Bn(u,v)  and 

B0  0 


|u||  =  (B  (u,u) 

B0  0 


1/2 


are  equivalent  to  (u,v)  and  ||uii  ,  respec 
tively.  Regarding  D,  we  suppose 


(2.3) 
and  that 

(2.4) 


0  <  D  (  u ,  u  )  ,  VO*  u  ^  H, 


I  u  ||  D  =  ( D  ( u ,  u ) 


1/2 


is  compact  with  respect  to  !!  *  !|  ,  i.e.,  it  has  the  property  that 

from  any  subsequence  which  is  bounded  in  !|  •  l|  ,  one  can  extract  a 
subsequence  which  is  Cauchy  in  II  •  II  D  . 

We  then  consider  the  variational ly  formulated,  self-adjoint 
eigenvalue  problem 


(2.5) 


Seek  '  (real)  and  0  *  u  H  such  that 
Bq(u,v)  =  vD(u,v),  V  v  H. 


Under  the  assumptions  we  have  made,  there  is  a  sequence  of  eigen¬ 
values 


and  corresponding  eigenvectors 


Ui ' U2 ' ' ' ’ 


which  can  be  chosen  to  satisfy 

(2.6) 


B  (  u  .  ,  u  . )  =  i  .  D  ( u  .  ,  u  . 

0  1  j  i  i  j 


=  o  . 

i  - z 


where  .  .  =  1 

i .  J 


for 


i  =  i  and  6 .  .  =  0  for  i  *  j 
J  1 3 


Fur  the: 


more,  any  u 


2.7; 


H  can  be  written  as 

V 

u  =  /  a  .u  . , 

j  j 

3*1 


with  a..  =  B 0 ( u ,  u ^. )  , 


where  (2.7)  converges  in 


The  eigenvalues  \  ..  satisfy 


following  well-known  variational  principles: 


Bq(u,u)  B(u1,u1) 

'  1  =  D(u,u)  =  D(ur'u3T 


(2.8) 


k 


min 

u-:H 


B  .  (  u ,  u )  B  ( u,  ,  u,  ; 
0  k  k 


D ( u , u )  D(uk,uk) 


k  =  2,3,.. 


Bo(u,ui)=0, 


i=l , 2 . , , , k-1 

(the  minimum  principle; 


and 


\ k  =  min 


max 


BQ (u,u) 


V,  --H  u-:V 


D (u , u) 


max 


bq(u,u: 


d  ( u ,  u : 


,  k  =  : 


k 


k 


dim  Vk=k 


(2.9) 


u--Uk=sp(u1 . uk 

the  minimum-maximum  principle). 


For  any  '  we  let 

(2.10)  M(  '  )={u:u  is  an  eigenvector  of  (2.5)  corresponding  to 

We  shall  be  interested  in  approximating  the  eigenpairs  of 
(2.5)  by  finite  element,  or,  more  generally,  Galerkin  methods. 


Toward  this  end  we  suppose  we  are  given  a  (one  parametric)  family 


{ S ^  }  o < y,  1  of  subspaces  Sv  H  and  we  consider  the  eigenvalue 


problem 

(2.11; 


Seek  ' 
50 


reai 


0  *  u. 


h  '  h 
2  (  u,  ,  v  )  ,  V  v  -c  Su. 


such  that 


The  eigenpairs  f • ^ , u,  !  of  (2.11}  are  then  viewed  as  approxima¬ 
tion  so  the  eigenpairs  (>  ,u)  of  (2.5).  (2.11)  is  called  the 

Galerkin  method  determined  by  the  subspace  for  the  approxima¬ 

tion  of  the  eigenvalues  and  eigenvectors  of  (2.5).  We  will  alsc 
sometimes  refer  to  problem  (2.11)  as  the  Galerkin  approximation  r 
the  problem  (2.5).  (2.11)  has  a  sequence  of  eigenvalues 


0  <  * 


h,  1 


h,  2 


Vn'  n  =  dim  Sh' 


and  corresponding  eigenvectors 


uh ,  1  '  Uh ,  2 . Uh  ,  N 


which  can  be  chosen  to  satisfy 


2.12 


. ,u 


h,  i  '  h, j 


=  \  . D  (  u.  .  ,  u, 

h  ,  i  h ,  l  h  ,  j 


i  -  J 


i  -  j  = 


. N 


Minimum  and  minimum-maximum  principles  analogous  to  (2.8)  and 
(2,9)  hold  for  the  problem  (2.11);  they  are  obtained  by  replacing 
H  by  and  letting  k  =  1.....N.  We  will  refer  to  them  by 

(2.8*)  and  (2.9*),  respectively.  Using  (2.8)  and  (2.9)  together 
with  (2.8*)  and  (2.9*)  we  see  immediately  that 

(2.13)  'k  'h  k'  k  =  I-2 . N  =  dim  Sh. 

For  every  k  we 


M ( '  )  =  { u : u  is  an  eigenvector  of  (2.11)  corresponding 

(2.10*)  n '  ^ 

t0  'h.k>- 

Because  of  (2.1) -(2. 4),  0  Is  an  eigenvalue  of  neither  (2.5) 

nor  (2.11).  It  will  be  convenient,  however,  to  introduce  the 

notation  * _  =  ' .  n  =  ° . 

0  h ,  0 

In  what  will  follow  we  shall  assume  that  the  family  { S,  : 

n 

satisfies  the  approximability  assumption 


2.14 


inf  ||u-f 


+0  as  h - >0 ,  for  each  u  ■-  H. 


It  follows  from  the  variational  principles  (2.8),  (2.9),  (2.3*)  , 


and  (2.9*),  and  assumption  (2.14)  that  ^ - »a ^ 


as  h - >0 


each  k. 


Our  analysis  employs  two  functions  $  (  '  )  and  •(>.(')  of  the 
non-negative  real  variable  \  which  are  associated  with  the 
eigenvalue  of  (2.5)  and  (2.11),  respectively.  We  define 


(2.i5: 


j*l -  2 - 


(2.16) 


^  ( \ )  =  min  ! 1  -  - - 

j  =  1  ,  .  .  .  ,  N  h  ,  j 


It  is  immediate  that  the  functions  are  non-negative  and  continuous 
in  '  and  that 

*(\)  -  0  if  and  only  if  <  =  \^  for  some  ; 


*  ( '  )  =0  if  and  only  if  '  =  ' ,  .  for  some  i . 

h  *  h, j  J 

In  the  following  lemmas  we  give  characterization  of  *  and 


-  10  - 


respec- 


$ 

tv 


w 


\ 


ft 

i 


r-* ; 

$ 


.Vj 


♦  ^  which  do  not  involve  the  eigenvalues  i  and  * 


h,  j* 


tively.  For  0  -  a  <  v  and  u,v  ■£  H  define 


(2.17) 


B ( ' , u , v )  =  B  (u,v)  -  'D(u,v] 


We  now  have 


Lemma  2.1.  For  all  0  .  i  <  i, 


2  .  18 


*  (  v 


inf 

u-:H 


sup  i B ( \ , u, v) 
v^H 


Hull  =1  Hull  =1 

Bo  Bo 


Suppose  A.  has  multiplicity  q,  i.e.,  suppose 


'k-1  *  'k  'k+1  lk+(q-l)  <  'k+q- 


Then,  for 

(2.19) 


n , 1  ^  1 

k-1  k 


k 


A  <  \ 


_  1  _  1  __  *1 
f  A  .  ‘  +  A  .  "  1  * 
i  k _ k+q  ! 


k+q  ^  2 


we  have 

(2.20) 


*  ( V )  =  11  -  —l 


k 


and 


'MM  =  sup  |  B  (  M  u ,  v ) 
v-:H 


(2.21 


II VII  =1 

Bo 


=  I  B  (  i  ,  u  ,  u )  I  ,  V  u  :  M(\jJ  with  !  u .! . 


Proo f .  For  u , v  : 


H  write 

\ 
u 

j  =  l 


=  >a.u.  ,v  =  >  b  .  v  . 

:  J  j  ; 

j=l 


Then 


(2.22) 


B ( ' ,u, v) 


■I 


a  .b  .  (  1 
1  J 


1  =  1 


2 .23 


sup  i  B  (  '  ,  u ,  v )  I  =  /  a.(l  -  - — ) 

vs  H  Lf“  J  vj  J 


II  v  ||  =1 

B0 


from  which  we  get 


sup  | B ( V , u , v ) |  =  inf 

vs  H  j  =  l , 2 ,  .  . 


1  -  7—  I  =  * (* 


II  U II  b  =1  I|V||  =1 

Bo  Bo 


This  is  (2.18).  (2.20)  follows  from  the  definition  of  f- (  '  ) 

an  examination  of  the  graphs  of  |1  -  |  for  j  =  1,2,... 

j 

(2.21)  follows  from  (2.20),  (2.22),  and  (2.23). 

In  a  similar  way  we  have 


Lemma  2.2.  With  H  replaced  by  S^,  by 


and  u. 


u^  ^ ,  Lemma  2.1  holds  for  $h(X).  (Relationships  analogous 
those  of  Lemma  1  will  be  indicated  by  an  asterisk.) 

If  is  an  eigenvalue  of  multiplicity  g,  then  *h 
'h  k+q-1  could  be  multiple  or  simple.  The  graphs  of  i>  (  '  ) 
#^('  )  is  given  in  Figure  2.1 


-l  +  X/Xk- 


■-1  +  X/X I 


Figure  2.1.  The  graphs  of 


v-  »*?.'  »*«v  iTi 


s‘,»7 


v. » . 


yi 


We  end  the  section  with  a  lemma  that  expresses  a  fundamental 


property  of  eigenvalue  and  eigenvector  approximation. 

Lemma  2.3.  Suppose  (\,u)  is  an  eigenpair  of  (2.5)  with  ||u||„  = 

'  D 

1,  suppose  w  is  any  vector  in  H  with  ltwllD  =  1/  and  let  \  = 
BQ(w,w) .  Then 


(2.24) 


r  -  \  =  II  w-u|l2  -  Mlw-Ullj. 

Bo  D 


(Note  that  we  have  assumed  u  and  w  are  normalized  with  respect 

to  ||  •  ||  here,  whereas  in  (2.6)  and  (2.12)  we  assumed  u^  and 

u.  are  normalized  with  respect  to  ||*||Q  .  ) 
h,  l  Bq 

Proof .  By  an  easy  calculation, 

II W-U||  2  -  Mlw-ull’J  =  l|w||2  -  2Bn(w,u)  +  null  2 


Then,  since 


-  VllwllJ  +  2  v  D  ( w ,  u )  -  v  Hull  ^ 


!l  w  ||  D  =  II  u  ||  D  =  1, 


iwiir  =  \, 
Bo 

I  UllJ  =  V, 

Bo 


BQ(w,u)  =  V  D ( w, u] 


we  get  the  desired  result 
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The  main  result. 


f  H 


1 


m 


For  i  =  1,2,...  suppose  is  an  eigenvalue  of  (2.5)  of 

ki 

multiplicity  q^,  i.e.,  suppose 

k.-l  k.  k.+l  k.+q.-l  k . +q .  k.+l 

111  2  1  2  2  2 

Here  k^  =  1,  k^  is  the  lowest  index  of  the  2nd  distinct  eigen¬ 
value,  k  is  the  lowest  index  of  the  3rd  distinct  eigenvalue, 

o 

etc.  Let 


(3.1 


r,  ,(h] 


inf  inf  ||u-r  Ils  ,  j  =  1 ,  .  .  .  ,q, 

u--M(V.  )  t-S.  c0 


Hull  =1 
Bo 

Bo(u'uh,k.)l=  •••  =  Bo{u'u 


„>  =  0 


The  restrictions  B  (u,u 


h,ki+j-2 


o 1 u ' uh . k . 1  =  ■■■  =  Bo|u’uh.k(+j-2 


sidered  vacuous 


if  q .  =  1.  Note  that  t  .  = 

1  1/1 


=0  is  con- 


.  ,  a  .  = 


t  ,  where  r  ^  and  are  the  quantities  introduced  in  Sec- 

i  i  i 

tion  1.  It  is  the  purpose  of  this  section  to  estimate  the  eigen¬ 
value  and  eigenvector  errors  for  the  Galerkin  method  (2.11)  in 

terms  of  the  approximabi  1  i ty  quantities  t  ..(h). 

1  <  J 

Theorem  3.1.  There  are  constants  C  and  hQ  such  that 

(3.2)  v,  .  ,  ,  C  y2  .(h),  VO  <  h  ..  hn ,  j  =  1 . q.. 

ki  +  j-l  h,ki  +  j-l  i ,  j  0'  -  - 


!  n  -3 
v  -  •  - 


'U  -U  1 

’  k.+j-l  h,k.+i-l"B_ 
2  J  2  J  0 


C';  .  .  ( h )  ,  V  o  <  h  h  ,  j  =  1 . q 

l  ,  J  0  : 


i  =  1,2,...  .  To  be  slightly  more  precise,  the  eigenvectors 


W  • 
(2.6)  )  . 


of  (2.5)  can  be  chosen  so  that  (3.3)  holds  (as  well  as 
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Proof .  Overview  of  the  Proof.  The  complete  details  of  the  proof, 


which  proceeds  by  induction,  are  given  below.  Here  we  provide  an 
overview.  In  Step  A  we  give  the  proof  for  i  =  1.  The  proof  is 
very  simple  in  this  case  and  rests  entirely  on  the  minimum  prin¬ 
ciple  (2.8*)  and  Lemma  2.3. 

The  central  part  of  the  proof  is  given  in  Step  B.  There  we 
prove  the  theorem  for  i  =  2,  proving  first  the  eigenvalue  esti¬ 
mate  (3.2)  and  then  the  eigenvector  estimate  (3.3).  In  particular, 
in  Steps  B.l  and  B.2,  estimates  (3.2)  and  (3.3),  respectively,  are 
proved  for  j  =  1 .  We  further  note  that  the  argument  used  in  Step 
B  proves  the  main  inductive  step  in  our  proof,  yielding  the  result 
for  i  =  i.  +  1  on  the  assumption  that  it  is  true  for  i  .  i.  Tc 
be  somewhat  more  specific,  in  Step  B.l  we  prove  (3.2)  directly  for 
any  i  a  2  and  in  B.2  we  prove  ( 3 . 3  >  for  i  =  _i  +  1  under  the 

assumption  that  ||  u,  .  -  u  . !!  _  - >0  as  h - ► 0  for  i  -  i  (cf. 

n  <  J  3  “o  ~ 

3.36)  . 


Details  of  the  Proof.  Throughout  the  proof  we  use  the  fact  that 

f  .  .(h)  can  also  be  expressed  as 

^  >  3 


>  ,  .  (  h )  = 

i  /  3 


3 . 1' 


inf 

u--M( 

ki 

(lull  =1 
B0 


inf  i| u  —  t  ,|n 
'-sh° 

V  ,uh,k. )=- • -=5C! 1 
l 


B0 (u,uh,k. 5 

l 


Step  A.  Here  we  prove  the  theorem  for  i  =  1. 

Step  A . 1 .  Suppose  (k,  =  l)  is  an  eigenvalue  of  (2 

K1  1 

multiplicity  q. ,  i.e.,  suppose 

1  12  q , 


5 )  with 


In  this  step  we  estimate  ' 


the  error  between 


-  \ 


„v 


S3 


h,  1  1  ' 

the  approximate  eigenvalue  among  ' 


h, 1 '  h,q. 


■  ,  ana 

that  is  closest 


to  \  ,  i.e.,  we  prove  (3.2)  for  i  =  j  =  l.  Note  that 


1 , 1 


(h)  =  inf  inf  Hu- »  !| 

u-.Mtij)  t;Sh  B0 

iiuiir  =i 

Bo 

is  the  error  in  the  approximation  by  elements  of  S,  of  the  most 

«  A 

easily  approximated  eigenvector  associated  with  \ 1 . 

From  the  definitions  of  ^ i  x ( ^ )  we  see  that  there  is  a  uv 

=  1  and  an  s,  a  S,  such  that 


a  M(( 

(3.4) 

Let 


with  ||  u 


h  B, 


i|uh  -  sh" B. 


h  ~  “h 


u. 


u  = 


■  sh  * 


By  the  minimum  principle  (2.8*)  we  have 


(3.5 


'h,i  "  '  l 


Bolsh,sh>  “ 


Now  apply  Lemma  2.3  with  (*,u)  =  ( *  , uv) ,  w  =  s.  ,  and 

1  *  *  il 

B0^sh,sh'’  This  yields 


(3.6) 


(3.4 


B0(sh'sh)  "  ' 1 


|lsh,uh''B0  'llish'uh,; 


VVb. 


O  ;  3  ,  U  ;  . 

h  h  so 


3.6)  yield  the  desired  result. 


Step  A . 2 .  In  this  step  we  prove  (3.3)  for  i  =  j  =  1.  Let 
be  eigenvectors  of  (2.5)  satisfying  (2.6).  Write 


ui ,u2 ' • 
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From  (3.6)  and  (3.7)  we  have 


'q  +1  •  ^ 

ql  j=qx  +  l 


I  !•;■’)■  ■  ■  I  (*!' 


j=q1+i 


=  lB(Wi'Vin 

*  1  81  lh,  l'“h,  1  ,uh,  l1*1 '  l*'h,  l)D<1;h,  :  ,uh,  :  '  1 

-  "h,l  *  Vh!l 


-  Cc‘(1(h) 


Hence 


luh,i  -  Za5‘’u; 


I  Hl 


)  1  1  o 
ill  c 


j-qj+i 


C(1  - 


"  a51)u 

z-  cL  .  U  . 

i  =  l  J  J 

Redefining  u.  to  be  — - - ,  we  easily  see  that 

j=l  0  J  Bo 

!! u  \\u  =  1,  so  that  (2.6)  still  holds,  and 

1  B0 


3.10) 


!!uh,i  -  Vb0  ■- 


as  desired.  Note  that  u  may  depend  on  h. 

Step  A . 3 .  Suppose  q  -  2.  From  (3.1)  we  see  that 
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ueM(^)  l-:Sh 

'lU|lB0=1  B0()'Uh,l)=° 

B0(u,uh, 1 )=0 


Choose  uh  €  M(^)  with  HuhHB  =  1  -  Bo^uh'Uh  =  0  and  Sh 
Sh  with  B0(sh,uh  2)  =  0  so  that 


(  3  .  12 


and  let 


ll'VShll80  =  f1.2<h)- 


uh  = 


'  sh  * 


/D(iV5h] 


'Dlvsh> 


Since  Bq(s^,u^  1)  =  0,  from  the  minimum  principle  (2.8*).  Lemma 
2.3,  and  (3.12),  we  have 


3.13) 


Kh, 2  "  X2 


-  V'b0  S  Ct?,2(h» 


This  is  (3.2)  for  i  =  1  and  j  =  2. 


Step  A . 4 .  In  Step  A. 2  we  redefined  .  Now  redefine  u2 , 

sc  that  u,  ,  . .  .  ,u  are  B_-or thogonal .  Write 
1  q1  0 


X(  2  ) 

a  .  u  . . 
1  3 

3  =  1 


Now,  proceeding  as  in  Step  A. 2  and  using  (3.13),  we  have 

11  -  vVi>  X  [a521]2  -  'IH2']2'1  -  '2’ j>' 


3=q, +i 


if 


1 


% 


r 


€ 

'm 


!« 

UH 


ft 

•to 

.j 

;gj 


'*.***, 

y> 


3.14) 


iUh,2  '  Laj  U  j "  B  Q  >  Cfl,2(h)- 


But  by  (3.10), 


al2’  =  B0(“h,2'“l 


=  B0(Uh.2'Ul-Uh,l) 


3.15 


5  l|Uh,2llB0l|UrUh,lllB0 


<-  C*  (h) 


1  C*l,2(h> 


Combining  (3.14)  and  (3.15)  we  get 


Z(  2 ) 

aj  U 

3  =  2 


j 11  Bq  s  11  uh 


,2  -  X^S'V  + 


5C£l,2(h)- 


{  2  ) 

a '  u 


i  =  2  J  3 

Redefining  u  to  be  -f -  ,  we  see  that  1|  u_  i|  _ 

gl  z 

Z  a(.2>uJ|R 

j=2  j  y  B0 


and  Bq(u1,u2)  =  0,  so  that  (2.6)  holds  and 


(3.16 


|uh,  2  "  U2!,B0  ■'  1  ,  2  (h)  ' 


which  is  (3.3)  for  i  =  1,  j  =  2 


Step  A . 5 .  Continuing  in  the  above  manner  we  obtain  the  proof  of 


(3.2)  and  (3.3)  for  i  =  1  and  j  =  l,...,q 


Step  B.  Here  we  prove  Theorem  3.1.  for  i  =  2. 
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/.  V.  /  \  y  V,  '*>  *  w*/  y. 


w  *•  fwwwmvi « 


Step  B . 1 .  Suppose 


q^+1)  is  an  eigenvalue  of  (2.5)  of 


V  <*2  * 


multiplicity  q  .  In  this  step  we  estimate  ,  the 

<-  ri ,  ic  2  ^  ^ 

error  between  \ ,  and  the  approximate  eigenvalue  among 

'  ,  ,  ,  that  is  closest  to  .  Note  that 

h,k2  h,k2+q2~l  k2 


3.17) 


2  ,  1 


1  (  h)  =  inf 


inf  !!u->  !! 


B, 


U-:M(h,  )  i-:S. 


l|u||R  =1 

B0 


Wrlte  'k2-i  =  ,lk2'  'h.k2-i  =  rh'h.k2-  and  'h,k2  =  Vi=2 


2w 


Then  0  <  r  <1  and  <p,  .  1  .  Let  i. . 

h  h 


h 


1+a 


-1 


.  From  (2.13)  and 


the  definitions  of  \.  and  v.  .  in  Lemmas  2.1  and  2.2,  respec- 

*■’  h ,  k 


tively,  we  see  that 
(3.18) 


•i  ,  <  V.  ,  . 

k2  h,k2 


A  simple  calculation  shows  that 


and 


>  —Ik. 

2  1+ir  2 


2<fih 

'h,k„  -  T  -T  ' k„  _  ^h'k  ‘ 
2  1+c .  2  2 


Hence  (3.18)  shows  that 


l+i 


Since  ' ,  . - •*  .  as  h - >0  (see  Section  2)  ,  we  see  that 

h ,  j  j 

2 

■>'  ,  - *1  ,  and  '  u - - - —  as  h - *0  .  Thus,  noting  thai 


n 


h  1  ~ 1 


l+r 


'—y  <  1  ,  we  see  that  we  can  choose  hn  such  that  0  <  h  <  h_, 


From  (2.19),  (2.20),  (2.19*),  and  (2.20*)  (see  also  Figure 

2.1)  we  get 


0  -  *h  ('h,kj  -  * 1 


=  i 


and  hence,  using  (3.19),  we  get 


' h, k2  ' k2 


(3.20 


h,k 

'h,k  | 

l  -  2  - 

l  2 : 

Vk 

2; 

' t 

K2  j 

v1  -  «,h1>’ 

r 

'  k  2  |_*  h  ( '  h ,  k  2  ) 

-  4>  M  \  |  ;  1 

■  h  ,  k  2  '  j  '  h 

- 

1  -  L4h('h,k2 

)  -  *  (  !  i  ’  1 

n'^2  J  “ 

'k2kirh.k2> 

-  1 

-  4>  f  \  v 

' ‘  h'k2  J 

Vrh, 


)  -  &  M  ) 

k  '  h ,  k„ 

C  £  • 


provided  that 
(3.21) 


1  2'k  *h(Vh,k  *  *  Mh,k 

L  £  £ 


[vrh,k  '  -  *'rh.k 'I'"1  -  1  2- 


We  will  now  show  that 


!3-22)  Vh.k,1  -  ♦"h.kJ>  •  Ci2,l(h;'  f°r  h  h0' 


where  C  depends  only  on  '  ,  ' 


,  and  '  ,  and,  in  particular 


3 

i 


is  independent  of  h,  and  h„  depends  only  on  >  ,  '  ,  > 

0  l  k2~:  .<2 

and  the  approximabi 1 i ty  of  the  eigenvectors  in  M('.)  by  Sv 

^2  h 

As  in  Step  A.l,  choose  u  ---  M('.  )  with  l|u,  ||_  =  1  and  s,  : 

n  x 2  n  Bq  h 

S,  such  that 
h 


13.23 


We  see  tha'1 


!luh  -  sh!lBQ  -  '‘2,:lh)  • 


is  the  B^-or thogonal  projection  of  u  onto 


S,  ,  i  .  e  .  ,  that 
h 


(3.24 


B0(uh  "  sh'V)  =  °'  V  v  “=  V 


Let  w,  =  u,  -  s,  .  Then 
h  h  h 


(3.25; 


(3.26; 


i|Wh!lBQ  =  '2,l(h) 


1  wh  !l  n  =  !IV!^ 

h  B0  h  Bq 


Next  we  write 


where 


sh  =  chuh  +  v 


B0luh-eh>  =  °' 


i.e.,  we  let  chuh  3c-orthogonai  projection  of  sx 


span  ( u^ ) .  Let  r^ 


1  -c,  )  u,  .  Then 
h  h 


which  implies 


(3.27; 


wh  =  rh  -  eh’  Bleh'rhi  *  °- 


i  I W  'I 

h  Bq  '  h'B0 


Furthermore,  using  (3.24),  (3.25),  and  (3.26)  we  have 
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1 


ch  =  VW 


'Sh»B0»"h»B0 


and 


c,  = 


'h 


Bo(sh'uh: 


=  B0(uh  "  wh'V 


=  »VS0  -  B0<wh 


1  -  !IW 


h  B, 


!  h ) 

2  ,  1  n  ' 


=  1 

>0,  for  h<h 


with  hQ  sufficiently  small.  Thus  we  can  assume 


0  <  c 


h 


Also  , 
(3.23 


To  see  this  we  refer  to  Figure  3.1  and  note  that 


il  r.  ||  _  =  !|  w  !l  cos 

h  Bo  h  B0 


and 


I w  ||  =  llu,[lQ  sin  i  =  cos  o. 

h  B0  h  Bo 


From  (3.28)  an  the  definition  of  r,  we  have 

h 


(3.29 


c.  1  II  w.  Iin 


and  hence 


(3.30) 


sup  IB( Vk  'sh'vl 


B0=1 


sup 

V---H 
!l  v  || 


B(‘h,k'sh'v 


=  [<=; 


Bo=1 

1  "  'h,k2'k2j 


I 

j  =  i 
*k. 


^I1 


h,k2  j 


2 


Combining  (3.27),  (3.29),  and  (3.30),  we  get 


sup 

v^S 


B('h,k2-Sh'v)l 


2  2  -  - 1  ? 

1  -  "Vb0>  -  'h,k2'k2!" 


ilvll 


h 

B0  =  1 


J-tT...'1  *  *  A'i 

«k„ 


!1  -  "V'b/!1  -  rh.k2-k2'2 


sup  (1  -  'h  k  V 

j  =  l,2 -  n'k2  k2 

*  k„ 


’h  B_  • 


,  2  2,2,- 

1  -  "wh“B  »  *  ('h.k. 


+  sup 
j=l , 2 , . . . 
-k„ 


-1.2 

h.k.'kj  'W 

i-  i- 


V  ‘  P 


=  Q. 


Thus,  using  (3.25)  and  (3.26),  we  have 
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»v*' 

».y 

yjy 

‘•ft! 

i*V 


4 


>S!i' 

»SM 


m 


s& 


ftil 


1 


ftd 

:ftj 


I 


..w 

1 


sa 


*  ( * 

"  h  h,k 


2  '  "  Sh ''  B , 


*  r  1 1  -  a  w.  II  ? ) + 2  ( \ 


!,Vb0>*  (1h.k2) 


n  1  ^  2 

+  SUP  (1  ~  rh  k  Vk1)2i|Wh"B  (1  ‘  11  Wh 11 B  ''ir  " 
i  =  i.2 _  n,k  k  n  b  h  b  j 


3  =  1,2 _ 

-k2 


"  t*2(rh.k  >  +  f2.1(h)]1/2' 


and  hence 


3.31) 


.  (h) 


4.  (  'w  u  )  -  'MV..  )  ,  — 

h  h,.<  h,k  24(V  , 

v  h ,  k 


where 


R  *  j-iTST../1 "  •  **‘rvk2». 

*kj 


One  can  easily  show  that 


24(1, 


Is  bounded  Independent  of  h; 


in  fact,  using  (2.20)  and  (3.19),  one  can  show  that 


(3.32 


24  (  v 


h,k2>  2V1-ocw2fl<h)]  ^<1-0 


provided  f  2  ^ ( h )  <  1/2,  i.e.,  provided  h  <  hQ  for  hQ  suffi¬ 


ciently  small.  Hence,  combining  (3.31)  and  (3.32)  we  get 


(3.33) 


*  h  ^ v  h  k  '  “  *  ^  v  h  k  ^  1  - ~ —  f  2  1  ^  ^  )  ' 

n  n,K2  n , k2  <  ( i-t  )  * ' 1 


which  is  (3.22).  Combining  (3.20),  (3.21),  and  (3.33)  we  have 


(3.34 


‘h.kj  '  ‘k2  -  C"2.1,h)'  f0r  h  -  V 


where 
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(3.35) 


C 


& 

$ 


ft 


& 

§ 


vis 

vg 

4 

!.»g 


vi1 

:& 

* 

•?> 

:*A 


1,3 


SS 


L3 


3 


iff. 

»vw 


^  (W 


This  is  (3.2)  for  i  =  2,  j  = 


Comment  on  Inequality  (3.34).  C  in  (3.34)  clearly  depends  on 


v  ,  v  ,  and  v ,  ,  but  is  independent  of  h.  Note  that  if  we 

■i  ^  ry  Kg 


were  considering  a  family  of  problems  depending  on  a  parameter 


we  could  bound  C  =  C ( r )  above,  independent  of  r,  provided  \ 


'  j.  (  f  )  was  bounded  above,  '  ^  ^  (  r  )  was  bounded  away  from 


and  r  -  r  ( t)  was  bounded  away  from  0  and  1.  It  follows  from 


3.31)  and  (3.33),  that  (3.34)  is  valid  for  h  ^  h_ ,  where  hA 

U  0 


depends  only  on  ,  \ ^  ,  and  the  approximability  of  the 


eigenvectors  ,  j  =  k2 . k2+q2~l  ,  by  Sh .  For  a  family  of 


problems,  h  ( r )  could  be  bounded  away  from  0  if  v.,(r)  was 


bounded  away  from  0,  (r)  was  bounded  above,  £  (  r  )  was  bour.d- 

k2 


ed  away  from  0  and  1,  and  the  eigenvectors  u.  =  u.(r),  j  = 

j  J 


k2 . k2+q2~l,  could  be  approximated  by  S^,  uniformly  ir.  r. 


Step  B.2.  Suppose,  as  in  Step  B.l,  that  has  multiplicity 

k2 


q2<  We  have  shown  in  Step  A. 5  that  we  can  choose  the  eigenvectors 


u, , u„ , 


of  (2.5)  so  that  (2.6)  holds  and  so  that 


(3.36 


u,  .  -  u  . ;  _ 

h,  j  r'B, 


C<  2  ■  (h)  ,  j  =  1 . q1  -  k2-l  . 


Write 


3 . 37 


V  '  2J 

h,k0  2-3  j  V 

‘  j  =  l 


From  (3.37)  we  have 


-  27  - 


which,  together  with  (3.34),  yields 


V1 


3.38)  I  V  |a 

3  =  1 


(k2h  2 


{1  '  Vk  /'i)l  + 
Ko  J 


I  f 


j=Vq2 


(k2}l 2 


:  -  '  ..  '  ,  ) 
^2  - 


a  C^2,l(h) 


Note  that  the  first  term  inside  the  absolute  value  is  negative  and 
the  second  is  positive.  In  addition 


C1  S  11  ~  'k  Aj*  -  C2 '  Vj  =  k2'  k2  +  1 . k2+q2-1' 


with  C1<C2  positive  numbers.  Hence  from  (3.38)  we  obtain 


V1 


(3.39 


I 


r  <^2h 2 


3  =  1 


,  Dl£Jfl(h)  +  D2 


I 


:=k2+q2 


r  ( ko ) i 

La  ‘  J 


and 


(3.40) 


I  - 

j=k2+q2 


(koh  2 


k2-l 


-  Vu®  +  =3I 


3  =  1 


r  (k2  2 

1  a 
L 


Write 


3.41) 


X, 


uh,i  -  ui  *  2_bi.jur  1  =  1 . V1  -  v 


.1  =  1 


Then,  by  (3.36), 


i 

s 

4! 


§ 


■i 


yd 

$ 


4'J 


3 


1 


S3 


3 .42 


Zb  ?  .  =  !!  u.  .  -  u . !! 

i  ,  3  h, i  l 


C<2,i(h)'  1  =  1 . V1’ 


Next  we  wish  to  find  constants  ,  so  that 

x  K  2  «*• 


V1 


3.43) 


:°i-I  *  ai 

3  =  1 


/  i  ■“  i  • 


Using  (3.41),  these  equations  can  be  written  as 
k2"1 


k2‘1 


u .  ,  >  (a  .  u  .  +  a  .  >  b.  ,  u ,  )  =  u  .  +  )  b.  .<»  . 

1  3  3  3  3  >'■  (  J  i  Z_  3,1  3 


0-44) 


=  a.  ,  i  =  1 , . . . ,k„-: 


Since  (2.14)  implies  £  ^  ^  ( h ) - *0  as  h - >0  ,  from  (3.42)  we  see 

that  the  b.  .  are  small  for  h  -  h^ ,  with  h.  sufficiently 
small,  and  hence  the  system*  (3.44)  is  uniquely  solvable,  and, 
moreover,  there  is  a  constant  L,  depending  only  on  k0  ,  such 


that 


3.45) 


k2’1 

I  V 

-3  =  1 


k  - 1 

1/2  -  2  . ,  .-'12 
.  v  r  k2h2! 

-  L  >  i  a .  i 


V  r 

Z.  !ai 

Lj=i 


Now,  from  (3.36)  we  obtain 


!Bo(uh.k2-V’ 

iBoluh,k  ■'J3  -  uh.j' 


l|uh,k2»B0l|u3  -  Vj\ 

llu;  ■  uh,3«B0 
o  J  .(h),  j  =  l . k2-i 
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Letting 


(3.46) 


V1 


(h) 

*2 


■  I  ‘S. 


we  see  that 


(3.47; 


1c  —  1 

2  ,,  .  1 / 2 
v'  r  •k2)'i2l 


I  K T  •  cvh)- 


and  thus,  from  (3.45! 


3 . 48 


k  - 1 

2  n  1/2 


I“i! 

3-1 


'  LCp.  (h) 
K2 


CP.  (h)  . 
k2 


Now  let 


3 . 49 


V1 


i.  =  u 


h,k  ]L  "  jUh,  j 
3-1 


Then  ». 


Furthermore,  from  (3.43)  and  (3.44)  we  get 


(3.50) 


W' 


fO,  i  -  k2  -  1 

j  k  -l 

la^2’-  Y  <i.b.  . 

i  1  jL.  J  J  .  1 

3  =  1 


From  (3.48)  and  (3.49) 


(3.51) 


h  ,  k  2  1  d  q 


-h,k2"Bo 

-V1  -  1  '  2 


|X"5 

'3  =  i  J 


Co  (h) . 

K2 
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Using  (3.34),  (2.21*),  (3.50),  and  (3.51),  and  the  fact  that 


i 

I 


(h) - >0  as  h - >0  ,  we  get 


„  Vh,k,  'k 

O'  ~  ,  (  h )  .  - -± - ± 

,  j.  V  i  > 

h,  k2 


"  W 


3 . 52 


B('k2<uh,k2'  in^- 


.  c-f  y  , 


Z\  ^  *■>  >  r  \  *s  J  y — *  >  * 

*.  K  2  “I  ^ 


=k2+q2 


where  C'  >  0  and  is  independent  of  h.  Combining  (3.42 


3.4  . ,  and  (3.52)  we  obtain 


*=k2+q2 


(k2h2 


c  r  2 ,  i { h ! 


k  -1 

(k„;  JL 


XV  ^  p  /  r — 1  n 

la, 


'"k2+*2 


V1 


c  k ;  ,  ( h )  + 

L  *  '  1 


I"V  I 


{"2’ 

,  a ,  i  o 


=  k2+<^2 


r  2 

C-2,l(h)  ^ 


k,-l 

2  '  v  '  u  \ 

Zl  V — >  V  /  r, 

L.ji  j  £  ia.  ‘ 

1  =  1  ''<=k2^2 


(3.53) 


!  b  .  i 


2  ! 1  2l 


l<=k0  +  q, 

C 


2,  1  (h)+  *  1  ,k  -1  !h)  ,/k2~1  L!  2L  [ai 

2  l.i  =  1 


‘2  V  2 


W-k2+q2 


°  ['■  2  ,  1  '  h  *  +  f 


•V1 


f  r  (k  u  2',  1  '2 

\L  h  I  I 

2  U=i  *■  J 


I 

=k2+q2 


r  \  .10 

T  ( k  ;,2‘  "  n 

f3.53;  is  a  quadratic  inequality  in  V  ,  a  2  T  i  "! 

L  1  J  1  J 

*-*  =k  +q 
2  ^2 


wncss 


solution  yields 


i 

1  =k„+q. 


k2_  1 


'2  l2  0  2 

J  -  -1 


2  v~  r  \ k0 ) -  2 

l,k2-l(h>Z.  [ai  ‘  J  + 
i  =  l 


(3.54) 


Combining  (3.39)  and  (3.54)  we  get 


V1 


2.  lai  j  >  Dl"2,l(h)  +  D2C^,k2-l<h>I.  [ai  J  +  D2C^' 


V1 


(k„K2 


and  thus,  since  r ^  ^  _^(h)  is  small  for  h  small, 


V1 


3.55) 


i  (•; 


( k2 ) ^ 2 


j  sV2,i(h)- 


Next,  combining  (3.40)  and  (3.55),  we  get 


3.56 


I  K 


(k2h2  2 

£  +  V2,l(h>- 


^k2+q2 


Finally,  from  (3.37),  (3.55),  and  (3.56),  we  have 


w1 


u,  ,  - 

h,k„ 


I  h  2  Ve„  -IK 


I  y1fa'^,|2+  5 


0  L  j=l 


i  k: 

j  =k_,  -*-q 


C2  2,l!h)’ 


W1  Ik. 


.i  a,  u 

j=k  J  J 

Redefining  u.  to  be  ■■  - - ,  we  see  that 

k2  k2+q2-1  : k  ! 

i!  -  a,  ti  .•  ;l  n 

j=k2  J  J  B0 


so  that  (2.6)  holds,  and 


3  .  57 


:|uh,k2  -  Vb0 


This  is  (3.3)  for  i  =  2,  j  =  1 


.''iV.’Vo'V1 


y  y  y.  y.  y.  c-’.--.  y.y,  y»y.  y.  -  .  y.y.  y.  < 


I, 


i; 


s 

■  J 


Comment  on  Estimate  (3.57).  In  the  proof  of  (3.57)  we  used 
(3.36),  which  was  proved  in  Step  A.  A  careful  examinations  of  the 
proof  of  (3.57)  show  that  we  that  we  did  not  use  the  full  strength 
of  (3.36),  but  only  the  weaker  fact  that  i|u,  .  -  u.i|_  - >0  as 

h.J  J  BQ 

h - >0  for  j  <  1c ^  —  1  ■  (Cf.  the  Overview  of  the  Proof.) 

Step  3.3.  Suppose  qa  -  2.  In  Step  B.l  we  estimated  'u 

'  .  In  this  step  we  estimate  '  ,  .  -  b 

k2  F  h,k2+l  k2+l 

We  proceed  by  modifying  problems  (2.5)  and  (2.11)  by  res¬ 
tricting  them  to  the  space 


=  (U  -Z  H  :  B0(u,uh)k2)  =  c; 


Sh  =  {U  -=  Sh  :  B^U.U^^)  =  0), 

h,k 

respectively,  i.e.,  we  consider  the  problems  (2.5  ^ )  and 

h,k2  h.k„ 

(2.11  )  obtained  when  H  and  S^  are  replaced  by  H  ^  and 


in  (2.5)  and  (2.11),  respectively.  (2.11 


has  the 


same  eigenpairs  (\,  .,  u,  .)  as  does  (2.11)  except  that  the  pair 

n  '  J  n  '  J 

h,k2  h , k„ 

(  '  h  k  '  uh  ic  )  iminated .  (2.5  )  has  eigenpairs  (  ' 


u.  '  which  in  general  depend  on  h.  Nevertheless, 


(3.58 


=  *,  .  '  =  0 . q  -2 

k  +'  k  +*  2 

2  2 


i-e.,  1  ,  the  eigenvalue  under  consideration,  is  an  eigenvalue 

K2 

h ,  k„ 


of  multiplicity  q2~l  for  Problem  (2.5 
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Its  eigenspace  is 


M 


( V  =  {u :  M(v  :  Bo(u'uh,k2 


=  0}  . 


We  can  now  apply  the  argument  used  in  Step  B.l  to  problems 


h,k. 


h ,  k. 


2 . 5 


)  and  (2.11 


and  we  obtain  (cf.  (3.34 


59) 


\.k2  +  l  “  'k2  +  l  ;  C,2,2ihl'  f°r  h  <  h0 


K  fc 

“ '  Ko 

Since  u^  ^  depends  on  h,  the  problems  (2.5  )  and 

h,k2 

(2.11  j  depend  on  h.  It  follows  from  the  Comment  on  Inequal 
ity  (3.34)  with  r  =  h  that  we  can  apply  the  argument  in  Step 


B.l  obtaining  C  and  hQ  that  are  independent  of  h.  To  see 

h,k  h,k_ 

this,  note  that  =  'k  '  by  (3.58),  »  '  'l'  the 

h,k  h,k  h,k 

minimum  principle,  and  t  =  f/i.  - *\ ,  /»,  =  /■  as 

k2-1  k2  k2_ 1  k2 

h - >0 ,  by  the  minimum-maximum  principle  and  the  fact  that  uv  , 

h ,  k 

- *u,  (cf.  (3.57)),  and  hence  V ,  is  bounded  from  above, 

^2  *2 


h ,  k. 


h ,  k 


is  bounded  from  below,  and  > 


2 

is  bounded  away  from  C 


h ,  k. 


and  1.  Then  observe  that  the  eigenvectors  in  M  *(',.) 


:ar. 


approximated  by  S^,  uniformly  in  h.  Note  that  in  Step  B.l  we 


also  used  the  fact  that  » 


h.  j 


ii  .  as  h - >0  for  j  = 


and 


It  is  easy  to  see  that  the  corresponding  fact 


-ue  : 


the  present  context . 

Step  B . 4 .  Suppose  q2  -  2  as  in  Step  B.3.  Here  we  show  tha 


u.  can  be  chosen  so  that  II  u,  . 

k2+l  h,k2+l 


uk  +1'*B  C'  2  ,  2(h:  • 
z  J 


Ws 


know  that 


(3.60) 


*  >  j 


!’  U 


h-j 


Cf 


'  9  ,(h).  j  =  q.+l  =  k 


Icf.  (3.16),  (3.14),  and  (3.57)).  Assume  that  u, 

\r  4-1 
K2  * 

have  been  redefined  so  that  (2.6)  holds.  Write 


61  ; 


U 


h ,  k  2  + 1 


I 

3  =  i 


(k  +i) 

a  .  u  .  . 

J  ; 


If  we  apply  the  argument  used  in  Step  B.2  to  u,  , 

h'k2+1 

we  let  k2  be  replaced  by  k2  +  l  and  use  (3.59)  instead 
(3.34),  we  obtain 


W1 


(3.62)  ||  u 


h , k2+i 


i  •; 


(k  ) 


u  ||  -  Cr  (  h '  . 

J  B0  2,2 


But,  by  (3.60) , 


k. 


=  B0 


*  !  BOluh,lc2'ulc2 


:lu.  -  u.  ,  '!_ 

k2  h,k2  Bq 


C.-  „  ,  ( h ) 


>2.2(h) 


u 


h,  k. 


and  hence 


j=k2+l  -  J 

Redefining  u,  ,  to  be  — - - -  ,  we  see  that 

-1 

<-  c  Z  v  ^2 

il  z  a,  u 

j=k„+l  J  3  B0 

t. 

!iukn  +  l"B0  =  2'  B0(Uko  +  l'Uj)  =  °'  3  =  1 . k2'  50  that  ;2‘6) 

2  0  2 

holds,  and 


(3.63  i 


Ay  y  *  U.,  .  •  _ 

h , k^+i  k2  +  l  B' 


which  is  (3.3)  for  i  =  j  =  2. 


Step  B . 5 .  Continuing  in  this  manner  we  prove  (3.2;  and  '3 
i  =  2  and  j  =  1 , . . . , q  . 

Step  C.  Repeating  the  argument  in  B  we  get  (3.2)  and  (3 


i  =  3,4, 


This  completes  the  proof 


Remark  3,1.  It  is  possible  to  use  an  alternate  argument  in  St 
B.l  i f  we  introduce  the  so-called  Riesz  formulas  for  the  spent 
projections  associated  with  an  operator.  We  suppose  the  space 
and  the  bilinear  forms  BQ  and  D  have  been  complexified  in 


usual  manner.  Let  P  and  P, 


be  the  B, 


lcgonal  pro 


^-1 

tions  of  H  onto  M ( ‘  ,  )  and  >  M ( '  ,  ,  .),  the  dir 

i_  h ,  kn  +  i 

*  i  =  0 


the  eigenspaces  M(» 


h , k  +i ' ' 


. ^2 


q  -1,  respectively. 


:duce  next  the  operators  T ,  :  H - >H  defined  by 


T  f  H 


B  Q  (  T  f  ,  v  )  =  D  (  f  ,  v  )  ,  V  v  K 


|BQ(Thf  ,v)  =  D(f  ,v)  ,  V  v  -i  Sh 


It  follows  from  (2.1),  (2.2),  and  (2.3)  that  T  and  are 

defined  and  compact  on  H.  Furthermore 


(3.64 


T-T  )fH 

h;  "B, 


C  inf  !!  Tf  -  HI 


It  is  immediate  that  (),u)  is  an  eigenpair  of  (2.5)  if  and  only 
_  1 

if  (//  =  v  J',u)  is  an  eigenpair  of  T.  Likewise  (  ,  u,  )  is  an 

n  n 


eigenpair  of  (2.11)  if  and  only  if  (/ 1 ^  ; 


is  an  eigenpair 


of  T,  .  As  a  consequence  of  (2.14),  T, - »T  in  the  operator  norm 


associated  with 


.  Let  f  be  a  circle  in  the  complex  plane 


centered  at  /.t,  =  ,  enclosing  no  other  eigenvalues  of  T. 

K2  k2 

Then  for  h  sufficiently  small,  I'  will  contain  the  eigenvalues 

Vk,+i =  4h!k.+i'  1  =  ° . q2'1'  of  v  Aiso'  v  and 

P,  ,  which  are  referred  to  as  spectral  projections  associated 
'  2 

with  T  and  /.t  and  T,  and  f.i,  ,  .,  i  =  0,...,q  -1,  respec- 

'  2  h  h , k2  +  i  ^2  ^ 

tively,  can  be  written  as 


(3.65) 


3 . 66 


'2  2vi 


h  ,  \  2  TFT 


(z-T)  1 dz 


(  z-Tu, )  dz  . 


These  are  the  Riesz  formulas.  With  these  formulas  we  can  derive 
an  eigenvector  error  estimate  which  will  lead  to  the  eigenvalue 
estimate  (3.34)  . 
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Let  u  with  u || Q  =  1.  Then  v.  =  u  ■ 

q2-l 

y  M  (  V  h  j.  +^)  and  from  the  formulas  (3.65)  and  (3.66)  we  obtain 
1=0  2 

II  u  -  V  II  =  II  (  p  -  p  )  u  II 

h" B  “ '  i  ,  rh,»,  U" B_ 

0  k2  k2  0 

*  [I27I  i  IZ-Th)'1(T-Th,<2-T)'!u  dzS 


3.67) 


1  I  —  1  rj 

rfr  >  (z-T.  )  ( T-T,  )  — — -  dz|| 

i/r  l  j  h'  h'z  -j.i 


f 7F  [27r  rad<r^  sup  «(z-Th)  II  x  H!T-Th:u; 

(-“k2  *  rad'r>  *  'Jh.k2*q2-irl|l:T-ThIUl!E0 
Cll  (T-Th)ullB  . 


(3.64)  and  (3.67)  yield 


!  3  .  68 


C  inf  !|T  u- 1  || 


=  C  inf  H/iu-t 


C  inf  !|  u-  t  H . 


'his  is  an  eigenvector  estimate:  it  shows  that  starting  from  any 


M ( ' 2 )  with 


2  =  1  we  can  construct  a  \\  =  ( u 


M !  *  K  ,  )  that  is  close  to  u.  We  now  use  (3. 5  S',  to  prove 


13.34). 


3y  the  minimum  principle  (2.8*)  we  have 


■r  ■  »’i 


(3.69) 


Bq{v,v) 


inf 

v"Sh 

llvtlD=l 

B0(V'Uh,i)=0' 
i  =  i . k2-l 


q2_1 


Since  v,  = 
h 


rM,  ,  we  know  that  B_(v.  ,  u,  .)  =  0,  i  = 

'h'k2  1  0  h  h,;L 


Thus,  from  (3.69)  we  find 


^  h  ,  k^  ^2 


f  Vh  Vh  1 

B°  fe'  '  V 


Combining  this  with  Lemma  2.3  and  (3.68)  we  obtain 


h,k2  k2 


_ - _ it  2 

'Vd  '|U|ID  B0 


k2  !IVd 


|uj|D  "D 


,  C||vh  -  u||‘ 


-  C  inf  ||u-l 


Since  this  is  valid  for  any  u  -s  M(' 


with  |ju||p  =  1  ,  we  see 


that 


'  h  k  "  '  k  -  C  inf  inf  iiu_<  yR 

'  2  2  u-:M(h  )  t-:Su  0 

k2  h 

Hull-  =1 

Bo 


=  C-'-2,l(h)' 

which  is  (3.34) .  We  note  that  the  proof  given  here  rests  on  equa¬ 
tion  (2.24)  in  Lemma  2.3  and  employs  formulas  (3.65)  and  (3.66) 

to  construct  v,  =  v,  (u)  that  is  B. -orthogonal  to  u,_  .  ,  i  = 

h  h  0  h ,  l 

0 . k2“l,  and  satisfies 
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We  have  already  seen  that  the  eigenvalue  estimates  in  Steps  A.l, 


A.3(...  can  be  based  on  Lemma  2.3.  Proceeding  as  we  have  here,  we 
see  that  all  of  the  eigenvalue  estimates  (3.2)  can  be  based  on 
Lemma  2.3. 
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Numerical  Computations. 


4  . 

In  the  previous  sections  we  have  analyzed  the  errors  in  the 
Galerkin  approximation  of  an  eigenvalue  problem,  concentrating 
especially  on  the  case  of  multiple  eigenvalues.  In  this  section 
we  consider  a  finite  element-Galerkin  method  for  the  approximation 
of  a  model,  one-dimensional  problem  with  multiple  eigenvalues, 
presenting  numerical  results  and  their  analysis  in  terms  of  the 
results  of  Section  3. 


Consider  the  eigenvalue  problem 


1  ^ 

u  U!j 


=  I  <p  '  (x)u,  x  I  =  (  -H  ,n  )  , 


(4.1) 


u( -n )  =  u( n ) , 


|X 

j 


■  fr-  u']|n)' 


where 

,  .  —<t  ,  ,  1  +o 

<p[x)  =  n  |x|  sgn  x,  0  <  a  <  l. 

It  is  easy  to  check  that  the  eigenvalues  and  eigenfunctions  -:'e  as 
shown  in  Table  4.1. 

Table  4 . 1 

Eigenvalues  and  Eigenfunctions  of  the  Eigenvalue  Problem  (4.1) 


i  i  .  u . 

i  i 


0 

0 . 0 

a. 

1 

1 .0 

cos 

<J>  (x) 

2 

1 . 0 

sin 

<?  (x) 

3 

4 . 0 

cos 

2<?  ( X) 

4 

4 . 0 

sin 

2 <p  (x) 
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1,2 ,3,4,  for  <'  =  .4.  These  errors  are  plotted  in  Figure  4.1  in 

log-log  scale.  We  clearly  see  the  different  rates  of  convergence, 

specifically  seeing  the  rates  h2  and  h1+°  =  h-1-4  for  the 

errors  in  \.  . ,  for  i  =  1,3  and  i  =  2,4,  respectively,  as 
n ,  r 

suggested  by  (4.2)  and  (4.3).  It  should  be  noted  that  the  esti¬ 
mates  presented  in  Theorem  3.1  are  of  an  asymptotic  nature  in  that 
they  provide  information  only  for  small  h  (or  large  n),  i.e., 
for  h  (or  n)  in  the  asymptotic  range.  From  Figure  4.1  we  see 
that  for  it  =  .4  we  are  in  the  asymptotic  range  quite  quickly, 
say  for  n  -  16. 

We  computed  u^  1  and  ^ •  the  approximate  eigenfunctions 


corresponding  to  and  ' ^  ,  respectively,  normalized  by 

||*|ID  =  1.  The  results  of  Section  3  suggest  that  u^  1  should  be 

close  to  C  cos  <*>(x)  and  u.  „  close  to  C  sin  ¥>(x)  (cf.  (4.2) 

h ,  2 

and  (4.3)),  where  C  is  such  that  C  sin  <P  ( x)  and  C  cos  '/>(>:) 

-1/2 


are  normalized  by 


=  1 ,  i.e.,  C  =  n 


To  illustrate 


this  point  we  computed  and  i  =  1>2,  so  that 


4 . 4 


f. 


K(i)  = 


u.  h  -  cj1 1  cos  <p{x)  -  sin  <p(x)|ls 


lu.  .  -  C  '  cos  2«>  ( x ) 

i ,  h  1 


is  minimal .  We  would  expect  that 


4.5 


r(2)  p ( 4 )  r(l)  -(3 
C1  ’  G1  ’  C2  '  C2 


and 

(4.6 


(1)  _  r(2)  _  r ( 3 )  _  (4 
1  C2  ~  C1  2 


C  =  . 564189583  .  . 


Table  4.2  shows  some  of  the  results  for  <<  =  .4.  We  see  clear  1* 
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The  increase  in 


the  results  predicted  in  (4.5)  and  (4.6). 
c|4^,  ci1  1  ,  and  ci,3,  with  increasing  n  is  due  to 

12  c. 

value  solver  we  used.  Table  4.2  also  shows  that  K ( 1 ) 
K(3)  <  K(4),  as  we  would  expect. 

The  last  column  in  Table  4.2  and  Figure  4.1  show 
ratios 


'  2  : 

1 

the  eigen- 
<  K ( 2 )  and 


that  the 


' h, i+1  -  V i+1 

X  ,  .  -  A  . 

h ,  i  i 


i  =  1,3, 


increase  as  h - >0 .  This  shows  that  in  the  whole  h-range  we  con¬ 

sidered,  the  approximate  eigenvalues  converging  to  a  multiple 
eigenvalue  are  well  separated. 


Table  4.2 


Numerical  Solution  of  the  Eigenvalue  Problem  (4.1)  for  =  .4 


n 

i 

'h,  i 

1  1  - 

K(  i 

‘I11 

c(i' 

°2 

_ 

" 

h  ,  i  - 1  i  -1 

' ,  .  -  y  . 

h ,  i  i 

i 

1 . 0716754 

.  2704 

0 

. 5637791 

0 

-  .  1124891 

-16 

1  .  5562955 

8 

2 

1  .  1 1  15481 

.3423 

0 | - .4151973 

-13 

.  5636998 

0 

■ 

"5 

5 . 0394692 

.  1075 

-t-1 

. 5558919 

0 

. 1317809 

-12 

1.1943249 

4 

5 . 2414639 

.  1191 

+  1 

.  5022638 

-13 

. 5516234 

0 

1 

1 .0175850 

.  1329 

0 

. 5641633 

c 

.  1596754 

— 

-12 

2 . 0041570 

■  16 

2 

1 .0352431 

.1881 

0 

- . 8916589 

-12 

.  5641519 

c 

, 

3 

4 . 2691915 

.  5259 

0 

.5636643 

0 

.  1  124328 

-13 

1 . 2575063 

4 

— 

4 . 3385100 

.  5869 

0 

- . 2689727 

-12 

.  5637697 

0 

1  .  0043740 

.6018 

X 

.  5641879 

0 

.  6411454 

-11 

-  ■  - 

2 . 600  3837 

:  32 

2 

1 .0113741 

.  1067 

0 

.  1323421 

-10 

. 5641830 

0 

3 

4 . 0666055 

.  2589 

0 

. 5641561 

0 

.  1970954 

-10 

1 . 4067517 

f 

4 

4 . 0936974 

.  3067 

0 

- . 7375504 

-10 

.  5641613 

0 

i 

i 

1 

1 . 0010921 

.  3305 

-1 

. 5641895 

0 

. 7729760 

-9 

3 . 5 1 9000 1 

;  64 

2 

1 . 0038431 

.  6202 

-1 

.8670648 

-9 

. 5641883 

0 

i 

i 

3 

4 .0166006 

.1289 

0 

. 5641875 

0 

. 3641341 

-10 

1  .  6437659 

i 

» 

4 

— 

4.0272875 

.  1653 

0 

. 1415775 

-8 

. 5641858 

0 

i 

i 

1 . 0002729 

.1651 

-1 

. 5641895 

0 

.  4535626 

-7 

4.9215830 

■  128 

2 

1  .  0013431 

.  3665 

-1 

.  3251219 

-7 

. 5641893 

0 

1 

2 

4 . 0041468 

.  6440 

-1 

. 5641895 

0 

.  4409247 

_  7 

2  .  0107071 

. 

4 

4 . 0083380 

.9135 

-1 

-.9  05611 

-8 

.  5641890 

0 

' 

1 

| 

1 

1 . 0000682 

.8255 

-2 

.  5641896 

0 

. 8070959 

-5 

7 . 0542522 

1 

i  256 

2 

1 . 0004811 

.  2193 

-1 

.7269570 

-6 

. 5641895 

0 

! 

3 

4 . 0010365 

.  3217 

-1 

.  5641896 

0 

.  6435344 

-6 

2  .  5706  705 

! 

1 _ 

JL 

4 . 0026645 

.5162 

-1 

- . 2601000 

-6 

.  564  1895 

o 

We  next  consider  the  case  when  <»  =  .01.  Table  4.2  presents 
the  same  results  for  w  =  .01  as  Table  4.2  does  for  a  =  .4. 
Figure  4.2  shows  the  graph  of 


,  vh,i+l-'i+l  .  ,  . 

log  - i - : -  -  i  =  1-3 

h,  i  i 


as  a  function  of  the  number  of  intervals  n  in  a  semi-logar  ithrr.i 

scale.  The  computed  values  are  indicated  by  o's  and  '  s .  The 

graphs  are  formed  by  interpolation  {solid  lines)  and  extrapolatio 

(dotted  lines).  We  note  three  related  phenomena  that  did  not 

occur  with  a  =  .4.  For  small  n  the  approximate  eigenfunction 

-1/2 

associated  with  i.  ,  is  u,  ,  -  »  sin  </>(x),  in  contrast  to 

h ,  1  h  ,  1 

- 1  '  2  -  i  '  2 

u.  ,  -  n  '  cos  <p(x)  when  o  =  .4.  We  remark  that  n  *  cos 
h ,  1  . 

_l/2 

( x )  is  more  easily  approximated  by  than  is  n  sin  «>  '  x ) 

for  all  0  <  o  <  1.  This  anomaly  is  present  for  n  -  64  and  for 
n  -  128  we  get  results  which  are  in  agreement  with  the  (asymp¬ 
totic)  results  in  Section  3.  For  „  and  i.  .  we  have  to 

take  n  i  256  to  get  results  which  agree  with  the  asymptotic 
theory . 

For  u  =  .01  we  see  that  K(2)  <  K(l)  for  small  n(n  _  64! 

and  K(2)  >  K(l)  for  large  n  and  K(4)  <  K(3)  for  small 

n(n.l28)  and  K(4)  >  K ( 3 )  for  large  n.  Recall  that  K(2)  v 

K(l)  and  K(4)  >  K(3)  for  all.  n  when  a  =  .4. 


Finally  we  note  that  when 

<  j  = 

.01  the  ratio 

\ 

h,  i  +  1 

i  +  1 

Ah,i 

%  . 
l 

1  —  1,0, 

first  decreases  as  n  increases,  then  for  some  n  the  two  eigen 
value  errors  become  equal,  and  then  the  ratio  increases  again. 


This  is  in  contrast  to  the  case  for  <;  =  .4,  in  which  the  ratio 
increased  over  the  whole  range  of  n  values.  We  further  note 
that  the  value  n  for  which  the  eigenvalue  errors  are  equal  — 
n  -  70  for  i  =  1  and  n  -  16C  for  i  =  2  (see  Figure  4.2'  - 
marks  a  transition  in  each  of  these  situations  from 
sin  <p  (  x ) 


-1  2 


to  u 


:o  Vi 

-1/2 


- 1  '  2 

n  '  cos  4>(x)  and  u 


h,  3 


h.  1 

-12 

n  sin  2 •*)(:-:! 


h,  3 


cos  2<o(x),  from  K  (  2  )  <  K  (  1  )  and  '  4  :  <  V. ' 


to  K ( 2 )  >  K ( 1 )  and  K(4)  >  K(3),  and  from 
1,3,  decreasing  to  increasing. 


h ,  i  + : 


We  have  thus  seen  that  for  </  =  .4  the  numerical  results  are 
in  concert  with  the  (asymptotic)  results  in  Section  3  for  the  who  1 
range  of  n  considered,  while  for  a  =  .01  they  are  in  disagree¬ 
ment  for  small  n,  but  are  in  agreement  for  large  n.  We  now 
make  an  observation  that  further  illuminates  these  two  phases  of 
error  behavior  —  the  pre-asymptotic  and  the  asymptotic.  Toward 
this  end  we  note  that  if  ('  ,u  ),  with  !! u  |!_  =  1,  and 

11  ID 

('h  l,uh  1  '  '  II  u^  i^d  =  1  '  are  firs,t:  eigenpairs  of  (2.5; 

and  (2.11),  respectively,  then 


(4.7) 


* .  1  -  =  llu.  -u  ||  -  \  II  u  .  — u .  Ii . 

h  ,  1  1  h  ,  1  1  Bg  1  h  ,  1  I 


=  inf  [  li  i  -u 


1  "B, 


'  l1!  1  Ul!lD' 


<!iD=1 


If  '  is  a  multiple  eigenvalue,  then  the  u^  in  (4.7)  can  be 
any  corresponding  eigenvector  with  Hu^ll^  =  1.  (Note  that  we  are 
here  assuming  u^  and  u^  ^  have  ||  •  !| length  equal  1,  whereas 


l 

i 

I 


in  (2.6)  and  (2.12)  they  are  assumed  to  have  (Ml  -length  equal 


to  1 . )  The  first  inequality  in  (4.7)  follows  from  the  minimum 


principle  (2.8*)  and  has  already  been  stated  in  (2.12).  The  firs 
equality  in  (4.7)  follows  immediately  from  Lemma  2.3  with  (•■  ,u)  = 


! W ■  w  =  Vi 


and  l  ®0'uh, 1 ,uh, 1 1  'h,l'  If 


with  ;i  t  ,)  =  1,  then  from  the  minimum  principle  (2.8*), 


4 . 8 


'  h  ,  1  '  1 


B0 ( 1  '  1 


Again  from  Lemma  2.3,  this  time  with  with  (',u)  =  (v,,u, ),  w  =  < 

-l  j. 

and  \  =  Bq ( <  , t ) ,  we  have 
(4.9|  B0(,.,)  -  q  , 

The  second  equality  in  (4.7)  follows  from  (4.8)  and  (4.9).  It  is 
clear- from  the  above  discussion  that  u^  can  be  any  eigenvector 
corresponding  to  \  . 

From  (4.7)  we  have 


4.1°)  -  q  *  »>-Vb0  -  V-VD'  V 


»  --  S, 


wi  th 


If  (  is  II  •  l|  _  -close  to  u,  ,  to  be  more  precise,  if  ■  is 
B0  1 

taken  to  be  the  Bn-pro jection  of  u,  onto  S.  ,  then  the  second 

U  IP* 

term  as  the  right  side  of  (4.10)  is  negligible  with  respect  to  th 

first  term.  This  follows  from  the  compactness  assumption  made  or. 

i*lln  in  Section  2.  On  the  other  hand,  if  j|u.,-ii|_  is  not  smal 
J  B0 

' h  i  ~  ' i  maY  still  be  small  because  of  cancellation  between  the 

two  terms  on  the  right  side  of  (4.10).  Regarding  the  case  ’  = 

.01,  this  explains  why  for  h  large  (the  pre-asymptot ic  phase  . 

-1/2 

we  can  have  uh  x  "  sin  </>  (x)  and  K(l)  >  K(2),  and  yet 

have  ' h  1 ,  the  approximate  eigenvalue  associated  with  uv  ,, 

-  50  - 


closer  to  than  Is  2'  ^e  aPProx^ma1:e  eigenvalue  assoc; 

- 1  •  2 

ated  with  u,  „  n  cos  'p(x),  while  for  h  small  'the  -cv: 

z\ ,  ^ 

- 1  '  2 

totic  phase),  we  have  u^  i  ~  ”  cos  ^  (x)'  K(l)  <  K ( 2 ) ,  anc 

' ^  1  closer  to  i ^  than  is  v ^  ,  showing  that  the  eigenvalue 

2 

error,  is  governed  by  inf  I! ) -u  |l  . 

n ,  l  l  _  _  IB. 

(--sh  0 

This  situation  is  very  similar  to  the  situation  with  r. 
fixed  and  “  varying,  as  can  be  seen  from  Table  4.4  where  compi 
tations  for  the  case  ,n  =  4  are  shown.  We  see  that  the  charac¬ 
teristics  observed  in  Table  4.3  regarding  dependence  on  n  are 
present  in  Table  4.4  regarding  dependence  on  -/  .  Namely,  the 
abrupt  switch  in  the  values  of  C^1^,  the  abrupt  switch  f: 

K(2)  <  K  (  1  )  and  K  (  4  )  <  K  (  3  )  to  K  (  2  )  >  K  (  1  )  and  K(4)  >  K  3  ’ 
and  the  abrupt  switch  from  decreasing  to  increasing  ratio  cf 
errors  near  the  parameter  value  corresponding  to  .  ,  =  ■  „  . 

a  ,  ...  . .  . 

mention  this  situation  —  <<  varying  and  n  fixed  —  since  i  -  : 
easier  to  understand  in  terms  of  perturbation  theory  (cf.  Kate 
[5])  than  is  our  original  situation  —  n  varying  and  fixe- 
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!  i 

■  1  1 

1 . 0520268 ! 

.  2338 

0 

.  8181940 

-11 

.  5634386 

0  i 

1 .0171143 

:  8 

1  2  : 

1 . 0529172  ! 

.2268 

0 

.  5645965 

0 

-  .  2916448 

-1 1  ! 

!  3  i 

4 . 8576239 ! 

( 

.  9593 

r> 

-  .  9346720 

- 1  3 

. 5597529 

A 

1.0164293 

!  4  ! 

i  1 

4.8717141! 

.9615 

o 

.  5604533 

0 

.  1167277 

_  i 

1 1 1 

1 - 

1 .0128661  i 

.1223 

0 

. 8717399 

-10 

. 5635957 

0 

1 .0111639 

16 

A  1 

1 . 0130098 ! 

| 

.  1052 

0 

.  5647369 

0 

-  .  8480131 

-9 

:  3  I 

4 . 2088367  j 

.  4650 

0 

.  2507177 

-10 

. 5636658 

0 

1.0C 87030 

f 

i  i 

!  4 

4 . 2106542 [ 

.4577 

0 

. 5642694 

0 

-  .  3101833 

-  1  C 

1 1 ; 

1 . 0032139 ! 

.  7274 

_  1 

-  .9345813 

-9 

.  5636031 

1 . 0 068764 

32 

< 

2 

i  i 

1 . 00 3  2  3 60  ! 

.  3568 

- 1 

. 5647430 

0 

.  1273043 

_  ^ 

!  3i 

4 . 0515675 ; 

.  2384 

0 

.  3745461 

-9 

.5638178 

1.0057284 

1 

i  4  : 

4 . 0518629 ! 

.  2205 

0 

. 5644172 

0 

-  .41  15544 

-  3  ! 

1  1  i 

,  ) 

' — 

1 . 0008063 ' 

.5369 

-  1 

-5 

.  5636032 

A 

1  .  C  C  1  7  3  6  2 

6  4 

1 

i2i 

1 . 0008077 ‘ 

1 

.  3398 

—  2 

.2462939 

_  - 

! 3 1 

4 .0128623  j 

.  1343 

01  .2743681 

-7 i  .5638240 

A 

1  . 0  0  3  5  99  7 

1 

!  4 

4 . 0129086 ! 

.9792 

_  1 
a 

.5544235 

r\ 

.3196172 

-  8 

i  1  i 

-  ( 

1 . 0002018 i 

.4  196 

__  ■i 

.5647430 

0 

.3356056 

_  c, 

l.C 064420 

i 

:  12  8 

;  2  i 

i  i 

1 . 0002031 ' 

1 

.  4775 

_ 

.7414162 

-6 

'  .5636032 

A 

i 

t 

1 3 ! 

4 . 0032196 ( 

.9166 

-1 

i  .2379072 

-6 

!  .5638239 

A, 

-  a  A  ,  a  e  ^  a 

i 

i 

j  4  | 

4 . 0032230 | 

.9745 

-2 

.5644235 

0 

.  1197135 

-5 

| 

1 

i1! 

1 . 0000504  j 

.4372 

_  i 

i  .5647429 

c 

.1061527 

—  ~* 

1.0216254 

‘  256 

i  o  ! 

1  Z  1 

1 . 0000515 

.4614 

-  iv 

-  .  1553659 

-4 

.  563603  1 

i 

A 

1  3  i 

4 . 0008054 1 

.5011 

-1 

.5644234 

A 

!  -  .  2  1  2  3  2  7  3 

-4 

.  0  -  **  - 

i 

!  4| 

4 . 0008079 ' 

.  7741 

-1 

!  .1165012 

_  5 

•  .5638238 

; 

.5638238 


Table  4.4 


Numerical  Solution  of  the  Eigenvalue  Problem  (4.1)  for  r.  =  4 

and  for  Various  <' 
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Further  information  may  be  obtained  from  Professor  I.  Babuska,  Chairman, 
Laboratory  for  Numerical  Analysis,  Institute  for  Physical  Science  and 
Technology,  University  of  Maryland,  College  Park,  Maryland  207^2. 


